Kadanoff— Baym equations and non-Markovian Boltzmann equation in generalized 

T— matrix approximation 
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. this, Sec. [II is devoted to the application of our approach to the T-matrix approximation. In Sec. IV we derive a 
non-Markovian Boltzmann equation in binary collision approximation. 
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A recently developed method for incorporating initial binary correlations into the Kadanoff- 
Baym equations (KBE) is used to derive a generalized T-matrix approximation for the self-energies. 
It is shown that the T-matrix obtains additional contributions arising from initial correlations. Using 
these results and taking the time-diagonal limit of the KBE, a generalized quantum kinetic equa- 
tion in binary collision approximation is derived. This equation is a far-reaching generalization of 
Boltzmann-type kinetic equations: it selfconsistently includes memory effects (retardation, off-shell 
T-matrices) as well as many-particle effects (damping, in-medium T-matrices) and spin-statistics 
C-j \ effects (Pauli-blocking). 
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I. INTRODUCTION 



Nonequilibrium properties of many-particle systems have traditionally been described by kinetic equations of the 
Boltzmann type. Despite their fundamental character, these equations have well-known principal shortcomings, e.g. 
I | (i) the short-time behavior (t < T cor - the correlation time) cannot be described correctly, (ii) the kinetic or the 
quasiparticle energy is conserved instead of the total (sum of kinetic and potential) energy, (iii) no bound states 
are contained, and (iv) in the long-time limit, they yield the equilibrium distribution and thermodynamics of ideal 
particles. 

An important generalization are the well-known Kadanoff-Baym equations derived by Kadanoff and Baym H , and 
Keldysh Q. However, the original KBE contain no contribution from initial correlations. Therefore, the KBE are 
unable to describe the initial stage of the evolution (to < t < r cor ) and the influence of initial correlations which can 
be important for ultrafast relaxation processes. 

To include initial correlations into the KBE, various methods have been used, including analytical continuation of 



|. A convincing 
initial state and 



| the equilibrium KBE to real times and perturbation theory with initial correlations [p|JlO|, 

£S| ■ solution has been presented by Danielewicz who developed a perturbation theory for a genera 
t-H ' derived generalized KBE which take into account arbitrary initial correlations. Finally, a straightforward and very 
intuitive method which does not make use of perturbation theory but uses the equations of motion for the Green's 
functions instead, has been developed in [0J^]. While perturbative approaches are restricted to situations where the 
coupling is weak, our method is valid for arbitrary coupling strength. In particular, it allows to consider systems with 
strong coupling, such as Coulomb systems at low temperatures and/or high density (e.g. metals and dense plasmas) 
and nuclear matter, and to include bound states. In Sec. || we briefly recall the main ideas of our method. After 



II. INITIAL CORRELATIONS IN THE KADANOFF-BAYM EQUATIONS 

Starting point of our approach is the first equation of the Martin-Schwinger hierarchy pd[| , 



(S ac — U ac )Gcb = 0~ab ± ihV a d,ceG ce ,bd> (1) 

/ d fr 2 V 2 \ 

with 5 ac = ih— + 5 ac , (2) 

V ot a 2m a ) 

together with an initial condition for G ce ^d, 

G ce ,bd\t c =t c =t b =t d =t = G c b(to)G e d(to) ± G c d(to)G e b(to) + C ce ^d(to)- (3) 
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Summation/integration over repeated indices is implied. Here, C denotes initial binary correlations in the system, 
and U is an external potential. The self-energy is defined by 

S QC G c fc = ±ihV a d,ceGce.bd = ±ifi>V a d,ce ^G c bG e d ± j • (4) 

Considering Eq. ([|) in the limit t = t' — ► to, we get explicitly 

dt^ac{ta, t)G c b(t, t ) — ±ihV a d,ce {G c b(to)G e d(t ) ± G c d(to)G e b{to) + C ce: bd(to)} . (5) 



Since the time integration is performed along the Keldysh-Schwinger contour, only time-local contributions of £ sur- 
vive on the l.h.s. The last term on the r.h.s shows that there must exist, in addition to the Hartree-Fock contributions 
(first two terms), another time-local part, which is related to initial correlations. That means, the self-energy has the 
structure (£ denotes the self-energy in the adjoint equation) 

S afc = Sf fc F + S o6 + £ a6 , (6) 

^ab = ^ab F + ^ab + ^ab \ CO 

with the time-local terms (here, we give the time arguments explicitly) 

Z I J(t,t') = X I »(t,t )6(t -t'), (8) 
t^(t,t') = t^(to,t')S(t-t ). (9) 

The further steps aim at the determination of these initial correlation terms and are sketched here, for details, we 
refer to Refs. Inserting (|) into (§), one obtains a Dyson-Schwinger equation for t,t' > to, 

(S ac — U ac — ^ac)G c b = <>ab, (10) 

which can be cast into the form G~^G c b — S a b- Functional differentiation of this equation with respect to the external 
potential U yields a Bethe-Salpeter equation for 8G/5U . Performing the same steps for the adjoint equation to (|l|) 
as well, a solution for SG/SU, which incorporates initial binary correlations, is obtained, 



UKj ab _ n n n 

— — LfarfLfcb T Lx ae — — 

OUdc OUdc 

where C has the time structure 



Gfb ± GaeGcfC e f,ghGgbGhd, (11) 



C a b,cd(t a tb, t c t d ) = C a b,cd(to)8(t a - t )S(t b - t )S(t c - t )S(t d - t ). (12) 



III. GENERALIZED T-MATRIX APPROXIMATION 



In the previous section we have obtained a formal decoupling of the Martin-Schwinger hierarchy by introduction 
of the self-energy. Furthermore, our approach shows, that initial correlations can, in principle, be straightforwardly 
included into this quantity. The next step on the way to a quantum kinetic equation is to choose a suitable approx- 
imation for the self-energy. Among the standard schemes are the random phase approximation (RPA), describing 
dynamical screening, and the T-matrix (or binary collision) approximation. The determination of E in these schemes 
without inclusion of initial correlations is well-known. For example, the T-matrix approximation leads to a non- 
Markovian Boltzmann equation. In Ref. |l2| ], this equation has been derived within the density operator technique. 
The nonequilibrium Green's functions approach, however, opens the possibility to derive two-time quantum kinetic 
equations with their well-known advantages (e.g. they fully include the kinetic and spectral one-particle properties). 
One-time equations are obtained by taking the time-diagonal limit of the two-time equations in a much simpler way 
than within the density operator technique. 

In the following, we will use the nonequilibrium Green's functions theory to derive a generalization of the usual 
T-matrix approximation, which includes initial binary correlations. 

According to Eqs. (|^Jl^), the self-energy is determined by the functional equations jl3| 
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E,,;, — — \ ,..„•., . { ScbGed ± S e bG c d + GcfGegC fgfihGhd i <"',./- 
S a 6 = ±ift ^ SaeGcd ± 5 ac Gde + G dgC 'agJhG f c Ghe ± 




E„ 6 = ±ift!4d lC e ■( G c fG e gC fqhhGhd ± C 



±l?l < GdgCagJhGfcGhe ± -Jy—Gfc > V C( 



(13) 

(14) 

(15) 
(16) 



Notice especially that, due to the structure of the self-energy, the arguments of the functional derivative in the 
equations for E and E are the same in both cases, 



E + t* IN = E + T, IN = E c + E /Ar + E /Ar = E. 



We now introduce an effective two-particle potential 5 by 

5E ao _ <5E ao SG e j 



SG, 



ef 



5U c d> 



(17) 



(18) 



dUcd SG e f 5U c d 
and define a generalized T-matrix p^ j, 

Tabbed — ^ab,cd ± ifr^ae,cfG fgGhe%b,hd ± ^ae,cfG fgGh e Cgb,hd- (19) 

In terms of Feynman diagrams, Eq. ( |l9| ) reads (the shaded block denotes the initial correlation C) 



± 
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± 





Comparing Eq. (|T^) with the solution ( p"l"| ) for 6G/6U, one obtains the relation 

<5E ab 



SU, 



cd 



±iHGdeTae,bfGf c . 



(20) 



So we could identify T with the correlated part of the two-particle function without the bare initial correlation C. 
The equation for the self-energy now takes the form 



V a d,ce {ScbG e d ± S e bG c d + G c fG eg C f gfihGhd + i^G cfG egTf g,bhGhd} 



(21) 



Functional differentiation of this equation yields a relation for S, which depends on T and on the quantity 5& N /5G = 
±ih$. Inserting this relation into Eq. ([H]), and evaluating the functional derivative 5E /JV /<5G, one arrives at two 
coupled equations for T and <!>, where self-energies and S have been eliminated. Keeping only the ladder-type terms, 
these equations can be written as 



T a b,cd — V a b,cd + &ab,cd + V a b,efG eg G fhC g h,cd + ifoV a b,efG eg G fh,T~gh,cd, 
&ab,cd = C a b,efGegGfhV g h,cd + ifr®ab,efG eg G jhVgh,cd- 

Eqs. ( |22| , p3| ) can be solved easily (see Appendix A), yielding an explicit expression for T, 

T~ab,cd — T a b.cd + 7a6,e/G eg G fhC g h,ijGikG jlT^i^d 

+ Tab,efG egG fhC g h,cd + Cab,efG e gGfh 



(22) 
(23) 



(24) 
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or, in terms of Feynman diagrams, 



+ 



+ 



+ 



Here, T denotes the well-known "ladder T-matrix" which obeys 

T a b,cd — V a b,cd + ifrV a b,efG eg G fhT g h,cd, 



(25) 



T 




V 


+ 


V 




T 



The system (^,£3|) can be regarded as a generalization of the usual T-matrix equation pq), where Eq. (|2J) shows 
explicitly the corrections which are due to initial correlations. 

If we now insert Eq. ( p^ ) into the equation for the self-energy (pl|), we obtain £ in T-matrix (binary collision) 
approximation, 

E ac = ±ihT a b : cdGdb ± ifiT a b ye fG eg GfhCgh,cdGdb 

± (ih) 2 T a b 7e fG eg GfhC g h,ijGikGjiTki 7 cdGdb, (26) 




T 


+ 


T 






+ 


T 








T 



and analogously £, 



Eac = ±ifrT a b,cdGdb ± iflCab,efGegGfhTgh : cdGdb 

± (ih) 2 T a b y efG eg GfhCgh,ijGikGjiTki yC dGd 



(27) 




Comparing these results with the predicted structure of the self-energies, Eqs. the time- local contributions 

are identified as 



or, diagrammatically, 




Interestingly, the correlation part E of the self-energy contains an initial correlation contribution, too, 
E QC = ±ihT a b,cdGdb ± {iKj 1 ' T abte fG eg G fhC g h,ijGikGjiTki,cdGd 

and, again in diagrams, 



rdb, 



(30) 




T 


+ 


T 








T 



With Eqs. (p6| p0[) we have found a generalization of the T-matrix approximation. In addition to the usual ladder 
term, the self-energies contain explicitly contributions of initial correlations. 

All relations derived so far are valid on the Keldysh-Schwinger contour. In order to obtain the Kadanoff-Baym 
equations and kinetic equations for the Wigner function, it is now necessary to specify the position of the time 
arguments of Green's functions on the contour. Then we obtain from the Dyson equation (|l^) the well-known 
Kadanoff-Baym equations for the correlation functions g^ (in the following, small letters denote quantities on the 
physical time axis, and the time arguments will be shown explicitly), 



di{sac(t, t) - «rf/(t, t)} g<(t, f) = / di{a> c (t, t) - a< c (t, t)} g<(f, t') 



In 



(If 



gl(t^{st b (i,t')-^ b F (i,t')} 



4c(*,*){s c t(*, t')-g> b (*,*')}, 
di{g> c (t, t) - g<(t, t)} &%(t, t') 



to 

The self-energies read in T-matrix approximation 

a%{t,t') = ±iht% !cd (t,t')gf b (t',t) ±iht% N d (t,t>)gf b (ti,t) ± i ht I a ^ cd (t,t')gi(t ,t), 
al(t, f ) = ±iM% tCd (t, t')gf b (t', t) ± iht% N d (t, t')gf b (t', t) ± iht a » cd (t, t>)g«(t>, t ), 
ag F (t,t') = ±ih (v ab>cd ±v ab4c )g%(t,t')S(t-t'), 
where the initial correlation contributions are given by 

Ccd^t') = [ dU^ ef (t,^g^ gh (i,to)c ghiCd (to)6(to-t'), 



fIN 

L ab,cd 



(t,f) = J dtc ab , ef (t )g^ gf Ah,tyXcd(t^y(to - 1), 

while the greater/less and the retarded/advanced T-Matrices obey the equations 

*l,cd(*.0 = *fi / div ab , ef g^ gh (t,^t% >cd (i,t') + ih J dtv ab , ef gf Uh (t,t)tf h!Cd (t,t'), 
0*> = V ^s(t - 1<) + l hJ div abtef g*{f h (t, t)tfl A cd (t, t% 

where we introduced the abbreviations 



(31) 



(32) 

(33) 
(34) 
(35) 

(36) 
(37) 
(38) 

(39) 
(40) 
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Se/,ihM = 9* /A (t,t')gfl A (t,t% Q%, gh {t,t') = g%(t,t')g%(t,t'), (41) 
g%f h (tj) = ±0 [±(f - <')] {^(t/) - QZ ftBh {tJ)} ■ (42) 
A further important relation is the optical theorem, which follows from Eqs. ( |39| ) and (pQ), 

*|,od(*'*') = ^^*a6 )e /(^i)ej, ffh (f,?)^, cd (l,t')- (43) 

Equations (|3l|-|4^) represent the Kadanoff-Baym equations in the generalized binary collision approximation. Here, 
the T-matrix contains contributions which are due to initial binary correlations. These additional terms can be 
separated from the "usual" T-matrix, and, in particular, do not influence the structure of the Lippmann-Schwinger 
equation (p0[). 



IV. NON-MARKOVIAN BOLTZMANN EQUATION 



In the previous section, we presented a far-reaching generalization of the usual T-matrix approximation by incor- 
porating initial correlations. This way, the Kadanoff-Baym equations have become sufficiently general to describe the 
evolution of a many-particle system on arbitrary time scales, in particular on ultra-short times after an excitation. 
Their solutions, the two-time correlation functions, contain a tremendous amount of information on the statistical and 
dynamical properties of a strongly correlated many-particle system, fully including damping (lifetime) of the one and 
two-particle states |p^| . However, in many cases the information contained in the Wigner distribution is sufficient. 
Therefore, in the following, we will derive an equation for this function, i.e., a kinetic equation in a narrow sense. 

For this purpose, we consider the Kadanoff-Baym equations (|l 32) in the limit of equal times t — t' and subtract 



them from each other. The result is an equation for the distribution function which reads, in momentum representation 
(we consider a homogeneous system without external forces) , 

t 

d C 

—f(p,t)=± J di {a > (p,t,i)g < (p,i,t)-a < (p,t,i)g > (p,i,t) 
to 

+g<{p, t, i)a> (p, t, t) - g> (p, t, i)&< (p, t, t)} 
= I(p,t) + I IC (p,t). (44) 

This so-called time-diagonal equation is a very general representation of a kinetic equation. The r.h.s. describes the 
influence of collisions as well as initial correlations on the Wigner distribution and is, in principle, determined by the 
exact self-energy and the two-time correlation functions. 

In order to obtain a closed kinetic equation, two major tasks remain: (i) an approximation for the self-energies 
has to be chosen, and (ii) the reconstruction problem, i.e., the determination of g< as a functional of the Wigner 
distribution, has to be solved. The first task has already been dealt with in the previous section, with the result being 
the generalized T-matrix approximation, given by Eqs. (|3"3|-|43|) . Let us now consider the reconstruction problem. In 
order to obtain the functional relation g^ = g<[f], we use the generalized Kadanoff-Baym ansatz (GKBA) proposed 
by Lipavsky et al. |1| , 

g*(p,t,t') =±{g R (p,t,t')f*(p,t') - f^(p,t)g A (p,t,t')} , (45) 

with f < = f and / > = 1 ± /. For the products then follows 

Sil(M') - G&(tJ)F&(f)+F*(t)g&(t,1?), (46) 

where we used the abbreviations = /^(pi)/^(p2) and (/12 = <?(pi,P2)- From Eq. ( f46| ) follow relations between 
the functions Q R I A and Q R I A which were defined in Eqs. ( p] , p| ), 

Sf 2 (M0 = £f 2 (M')^i2(O, (47) 
(^(MO^-^W^M'), (48) 

where we introduced the Pauli blocking factor ~N\i = 1 ± /(pi) ± 7(p 2 )- 



G 



Now we insert the self-energies in T-matrix approximation, Eqs. (|33|-|38|), into the time diagonal equation ([44), 
replacing t< with the help of the optical theorem (|43| ) and by means of the reconstruction ansatz (|f6|). The result 
is the collision integral /, 

dp 2 dpi dp 2 



i( Pl ,t) = (my 



dt dt dt 



Fn(?)F<$) - F<(t)F>(I) 



{2ttH) 3 {2nh) 3 {2-Khf 

{t*(pip 2 t, PlP2 f)g R {t, = t)t A {pip 2 % PlP2 1)^(1,;) 
+t R (p 1 p 2 t,p 1 p 2 t)g^(ii)t A (p 1 p 2 ip 1 p 2 i)g^,t) [F>(t)F<(E) - F<(t)F>Ct) 

-5f 2 (t, t> fl (pip 2 t,Pip2 1)^ 2 (l, I)i A (piP2 1 pipat) [^©^iK?) - F<(i)F>(t) 
- g* (t, t)t R ( Pl p 2 i, p!p 2 f)(?f 2 (I, (l>lp2 l P1P2 1) [F& {t)F< (I) - F< (t)F> (I)] } 



(49) 



with g^ 2 A = g R ^ A (pi,p 2 ) and t' R ^ A (pip 2 t, Pip 2 t) = (pip 2 \t R l A (t, t)\ P1P2), and the collision integral arising from 
initial correlations I IC , 



ih 



dp 2 dpi dp 2 



dt 



(2ttH) 3 (2Trh) 3 (2nh) 3 

X {^(PlP2^,PlP2?)^(PlP2F,PlP2 - £(PlP2* ) PlP2*)fr A (PlP2*,PlP2*)} 



(ihf 



dp 2 dpi dp 2 dpi dp 2 



dt dtdt 



with 



(50) 
(51) 



(27rfr) 3 (2ttH) 3 (2tt/i) 3 (2tt/i) 3 (2irhf 
x {* fl (PiP 2 *, P1P2 *)£(piPa i, p x p 2 l)t A (p!p 2 1, pip 2 1)^(1)6^(1, t) 
+g? 2 (t,i)N 12 (i)t R {p 1 p 2 t,p 1 p 2 I)lC(p 1 p 2 1, ^(PxPal, P1P2 t)} 

£(piP2 *. PlP2 *) = io)c(piP2, P1P2; *o)^12(*0,*)- 

With Eqs. (HHH) 

we have obtained a very general quantum kinetic equation. The character of its approximations 
goes far beyond that of the usual Boltzmann equation. The collision integral /(pi,i) was derived without any 
approximation with respect to the times and thus fully includes retardation and memory effects which is usually 
referred to as non-Markovian behavior. Many-particle effects, as for instance self-energy and damping jl4j], and spin 
statistics effects (Pauli blocking) are included. So far, no restriction has been introduced with respect to the retarded 
and advanced propagators g R l A . In principle, they are to be determined self-consistently from their KBE which 
follow easily from Eq. (|l^). However, to avoid this essential complication, in most cases approximations are used. For 
example, in the quasiparticle approximation, the propagators are given explicitly by 

g? 2 /A {t, = -i_e[±(t - f')] e i[£i2+ir 12 ](t-t') (52) 



[ihf 



with E a = ^ + f£ + Re of + Re of and T 12 = Im of -r «u«, 2 • 



Ira a? 



Furthermore, the retarded and advanced T-Matrices are many-particle generalizations of the familiar T-Matrices 
of quantum scattering theory. They have to be determined from the Lippmann-Schwinger equation ( |40| ) which reads 
in momentum representation 

* flM (piP2 1, p[p' 2 = »(pi - p'JpnhfSfa + p 2 - pi - p' 2 )S(t - t') 



+ih 



dpi dp 2 



(27Tft) 3 (27T^) 3 

^(p 1 ,p 2 ;t,^(p 1 p 2 f,p' 1 p 2 t' 



dtv(px - pi)(2ttS,) 3 (5(pi +p 2 - pi - p 2 ) 



(53) 



The collision integral I IC (pi 1 t) contains the terms arising from binary correlations, existing in the system initially. 
It should be stressed explicitly that the structure of these contributions is completely general and does not depend on 
parameters characterizing the system, such as coupling strength or degree of degeneracy. Furthermore, the inclusion 
of initial correlations does not depend on their actual form, i.e. the form of the function c. The damping of the 
two-particle propagators leads to a decay of this collision term, i.e. the initial correlations die out on a time scale 
which is determined by the one-particle damping rates Jl4| . 

Finally, we want to remark here that our result for the non-Markovian Boltzmann equation is in agreement with 
the result derived within the framework of the density operator technique, see [121. 
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APPENDIX A: SOLUTION OF THE GENERALIZED T-MATRIX EQUATIONS 

We rewrite Eqs. (E2LE3J), which contain the ladder- type terms of the generalized T-matrix (19), 



%b,cd — V a b,cd + &ab,cd + V a b,efG eg G fhC g h : cd + i^V a b,efG e gG}hT g h,cd, 
&ab,cd = C a b,efG eg GfhVgh y cd + ifr®ab,efG eg G fhV g h,cd- 



Due to the structure of Eq. (Al), T can be split into three parts, 

* ab.cd 



ab.cd 



ab.cd 



r (c) 

ab.cd' 



n-{A) 
ab.cd 



V a 



b,cd 



ihV ab,efG eg G fh T g { ^ cd 



T ab B ld = Vab,efG eg Gf h C ghlCd + ihV a b,efG eg G fh.T g ^ cd , 

r ab,ld = ®ab,cd + ihV abief G eg Gf h T g { ^l d . 



(Al) 
(A2) 



(A3) 
(A4) 
(A5) 
(A6) 



Obviously, Eq. (A4) coincides with the well-known ladder equation of the T-matrix approximation. Thus, we can 
ident ify with the usual T-matrix T. The ladder equation (A4) now serves as a basis for the solution of (A5) 

and (A6). If one assumes for T' s ) the form 



n~(B) rp *~i *~i s~1 

1 ab,cd — 1 ab,ef^ eg^ fh<^ gh,cd, 



(A7) 



Eq. (A5) is valid if (A4) holds. In order to determine T^ c \ Eq. (A2) has to be considered. This equation is fulfilled 
if $ is of the structure 



&ab,cd — C a b,efGegGfhTgh t cd, 



(A8) 



if the adjo int equation to (A4) is valid. Due to the symmetry properties of T, Eq. (A4) and its adjoint are equivalent. 
Inserting (A8) into Eq. (A6) and assuming T^ G > to be of the structure 



1~ab cd ~ G a b,efG e gGfhTgh.cd + i^T a i, te fG e gGfhCgh,ijGikGjlTkl } cd, 



Eq. (|Aq) is fulfilled, again under the assumption (A4). Collecting all parts together, T can be represented as 



T~ab,cd — T a b yC d + ihT a b^ e fG eg G fhC g h,ijGikGjlTkl !C d 

+ T a b^efG e gG fhC g h,cd + G a b,efG eg G fhT g h : cd, 

together with the equation for the well-known "ladder T-matrix" , 

T a b,cd — V a b,cd + ifoV a b, e fG ' eg G fh,T g h,cd- 



(A9) 

(A10) 
(AH) 
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